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ABSTRACT 

Cooling of primordial gas plays a crucial role in the birth of the first structures in 
our Universe. Due to the low fractional abundance of molecular species at high red- 
shifts, spontaneous emission rather than collisions represents the most efficient way 
to cool the pristine plasma. In the present work, radiative cooling functions are eval- 
uated for the diatomic species HD, HD"*", HeH+, LiH and LiH+. Cooling functions 
for the triatomic ions Hj^ and H2D+ are also considered. Analytic fits as functions of 
temperature are provided. 
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1 INTRODUCTION 

In its present shape, our Universe appears to be a large 
ensemble of bound structures. These are the results of an 
intricate net of chemical and physical processes of increas- 
ing complexity, the first of w hich i s frag mentation of the 
primordial gas clouds, see e.g. ISilkl (| 19831 ). One of the fun- 
damental ingredient of this ensemble of processes is the cool- 
ing of pristine plasma. Indeed, the birth of the first objects 
during the evolution of the primordial gas could not have 
occured without significant cooling, as the gas temperatures 
involved in the initial phase of the Universe are too high 
to allow the early clouds to gravitationally collapse. Among 
the plasma constituents, molecules play the most significant 
role as coolers of the pristine gas. Molecular cooling takes 
place mainly via rotational-vibrational (ro- vibrational) tran- 
sitions between the internal, nuclear motion degrees of free- 
dom. For temperatures below about 8000 K molecular cool- 
ing processes are mu ch more efficient tha n the corresponding 
atomic ones, see e.g. iMiller et all (|2010h . 

In models, the standard way of capturing the effects of 
cooling by a given species is via the cooling function, which 
gives the energy lost per second at a specified temperature. 
Cooling functions are made up by two main contributions. 
First, molecules spontaneously emit radiation while making 
ro- vibrational transitions; second, collision induced emission 
can occur; this is particularly important for high symmetry 
species where spontaneous emission can be very weak. 

In a plasma, two key features render a molecule an effi- 
cient cooler at low temperature: its fractional abundance 
and its dipole moment. In the chemistry of the primor- 
dial Universe H2 represents the most abundant molecular 
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species. Its role as an efficient primordial coolant has been 
widely described, e.g. in the formation of the first stars (e.g. 
Abel et al. 2002) and galaxies (e.g. Bromm et al. 2009, Ben- 
son 2010). Both coUisional and radiative H2 cooling con- 
tributions have been considered in the literature. The role 
of different coUisional partners (H2, He, H^ and e~) was 
considered by Glover & Abel (2008). Together with H2, the 
cooling properties of its cation H J due both to radiative and 
coUisional (H and e~ impact excitation) pathways have been 
described (Suchkov & Shchekinov 1978). Both H2 and HJ 
are customarily included in molecular cooling models be- 
cause of their high relative fractional abundance. However, 
because of their lack of a permanent dipole moment, only 
weak electric-quadrupole or magnetic-dipole transitions are 
allowed in these species, severly limiting their cooling effi- 
ciency, particularly at low density. 



Considering the mole fractions of chemical species given 
by chemical networks such as the one by Galli & Palla (1998) 
(hereafter, GP98), other molecular species should be added 
to the list of coolants. Grouping the primordial molecules as 
deuterated or helium-containing or lithium enriched, HD, 
HD^, HeH"*", LiH and LiH"'" represent the candidates to fo- 
cus on. Their fractional abundances are expected to have 
the following freeze-out values {z ~ 10): /(HD) ~ 10 
/(HD+) ~ 10 
/(LiH)+ ~ 10- 
at high z, but the most likely species are H;^ and H2D+; 
their fractional abundances are thought to be approximately 
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It is hard to form triatomic molecules 
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lisional contributions to the cooling have bee n evaluated in 
difLfe rent density regimes for HD ad LiH jLipovka et al.l 
(j2005l ). GP98). In particular, the effect of the most abun- 
dant collision partner H on this cooling pathway has been 
investigated. However, in order to complete this description. 
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the contribution of radiative cooling has to be taken into ac- 
count as well. 

The present work is organized as follows. In Section [2] 
the equations used to compute the radiative cooling func- 
tion are presented, describing the methods adopted for each 
molecule. Some details on previous calculations on coUisional 
cooling functions are given. Section[3]summarizes our results 
and provides fits to the calculated data. 



2 RADIATIVE COOLING FUNCTION: 
DEFINITION AND CALCULATIONS 

The total cooling function for the chemical species X is 
usually approximated as (e.g. HoUenbach & McKee 1979, 
GP98): 



Ax 



1 + 



n(collider) 



(1) 



where Ax.lte represents the local thermodynamic equilib- 
rium radiative contribution to the cooling function, n'^'^ and 
n(collider) are the critical and collider density, respectively. 
The former is defined by: 
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Ax(n(collider)) 
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Ti(collider)— 



where Ax (n (collider) — >■ 0) is the low-density limit of the 
cooling function. Eq. ([!]) is valid for all density regimes of 
the collider particle. In the high-density limit Ax reduces 
to Ax.LTE, implying that in this regime local thermal equi- 
librium is reached and that cooling is dominated by ro- 
vibrational transitions of the molecular species X. 

Denoting with the letters "it" and repectively the 
upper and lower ro- vibrational states between which a spon- 
taneous transition may occur, the radiative cooling function 
is defined as: 



W{T) = 



Ei){2Ju - l)gun^ (3) 



where Z{T) represents the partition function of the 
molecule, Aui the Einstein coefficient for the u — >■ I tran- 
sition, Ju the rotational angular momentum of the upper 
rotational state, Qu the nuclear spin degeneracy (equal to 
1 for all calculations below), Eu and Ei the energies of the 
upper and lower state, and n„ the population of the upper 
level. As discussed above, under LTE conditions it holds: 



W{T) = ALTE(r) 

Z(r) = ;^(2J. + l)3,e- 



(4) 



Radiative cooling functions have been computed in the 
present work using the equations above. As kg^ — 
1.43878 cm K, the Boltzmann factor n„ implies that at a 
temperature of 1000 K the cooling function is mostly deter- 
mined by transitions with Eu < 700 cm~^. 

Data on the emission probabilities and energy levels of 
the molecular species under analysis are therefore required 
in order to compute Alte. Depending on the available data, 
different strategies have been adopted for each molecular 
system considered. For HD and HcH^, the relevant data 



have been calculated respectively by Abgrall et al. (1982) 
and Engel et al. (2005). To obtain the necessary data for the 
other diatomic molecules under investigation the program 
LEVEL 8.0 by Le Roy (2007) was used. This program solves 
the radial one-dimensional Schrodinger equation 

fi2 d2*(r) 



2fi dr2 



+ Vj{r)<f{r) = S„,j*(r) 



Vj{r) = V{r) + 



fi^[j(j + 1) - n'] 

2/ir2 



(5) 



where jj, is the reduced mass of the system, J the total ro- 
tational angular moment, Q the projection of the electronic 
angular momentum along the internuclear axis (57 = for 
the states of our interest) and V{r) is the diatomic potential 
energy curve, given as input. Once the ro- vibrational wave- 
functions *I'^,j have been obtained level can also calculate 
the Einstein coefficients by 

Av' J'^u" J" ~ 

167V'^\ S(J',J") _3 , , , ,,,,,,, , ,2 (6) 

where S{J' , J") are the Honl- London rotational intensity 
factors, resulting from integration over the rotational de- 
grees of freedom, v is the wavenumber of the transition, 
Af(r) is the diatomic dipole moment curve and is given as 
input. In the common case where i> is measured in cm~^ and 
Al in Debye the constant in brackets in Eq. (|6]) assumes the 
value 3.1361891 x 10"^ cm^ D"^ s"\ 

Calculations were carried out for the electronic ground 
state of the molecular systems under investigation. The fol- 
lowing subsections examine each molecule considered and 
give detailed information on the potential energy and dipole 
moment curves used to perform the calculations. Compar- 
isons with existing data are also provided where possible. 



2.1 HD 

The role of HD in cooling phenomena of primordial plasma 
has been widely examinated in several studies (Flower et 
al. 2000, Lipovka et al. 2005, Ripamonti 2007). The ef- 
fect of HD cooling on the formation of massive primordial 
stars was analysed by Yoshida et al. (2007) and McGreer & 
Bryan (2008). Despite its lower fractional abundance, HD is 
a more efficient coolant than molecular hydrogen H2, espe- 
cially for T < 2000 K (GP98). This is because HD molecules 
possess a small permanent dipole moment (~ 8.3 x lO^** D, 
Abgrall et al. 1982) due to the asymmetry of the nuclei; by 
contrast, in non-deuterated molecular hydrogen only very 
weak quadrupole transitions are allowed. 

Lipovka et al. (2005) provide a two-parameter fit of 
the coUisional cooling function, depending on temperature 
and H fractional abundance. Here we evaluate a radiative 
cooling function which takes into account both dipole and 
quadrupole electric transitions, for which we used the calcu- 
lations by Abgrall et al. (1982). Results are shown in Fig. [TJ 
and compared with the high-density limit HD cooling func- 
tion fit provided by Lipovka et al. (2005), where the contri- 
bution of the lowest ro-vibrational levels up to w = 3 were 
considered. A deviation from the high-density limit coUi- 
sional contribution to cooling function is found for tempera- 
tures above 700 K, where highly excited ro-vibrational states 
are expected to be populated. 
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Figure 1. HD radiative cooling function. Circles: dipole contri- 
bution, present calculation; blue dotted- dashed curve: quadrupolc 
contribution; red dashed curve: high-density limit fit by Lipovka 
et al. 2005; green solid curve: fit given in Table [T] 



2.2 HD+ 

The HD cation and its coUisional contribution to the cooling 
during the evolution of primeval plasma has been examined 
by Glover & Abel (2008), where calculations were carried 
out assuming that there were no significant differences be- 
tween the molecular ions and HD"*" in the high-density 
cooling limit. These authors provide an analytical expres- 
sion for Alte, calculated using molecular data from Karr & 
Hilico (2006) on HD+ (energy levels) and Peek et al. (f979) 
on Hj" (quadrupole moments). 

Our calculations of the HD"'" radiative cooling func- 
tion were carried out using the potential energy and dipole 
moment curves provided by Esry & Sadeghpour (f999), in 
which an adiabatic reformulation of the HD^ Hamiltonian 
that recovers the isotopic splitting of electronic states is pre- 
sented. The resulting cooling curve is shown in Fig. [2l and 
compared with the fit by Glover & Abel (2008). It should 
be noted that our calculations predict significantly increased 
cooling at temperatures below 100 K, where the cooling by 
HD+ rotational transitions is very much more efficient than 
the one due to its parent ion. 



2.3 HeH+ 

HeH"*" is the first molecule formed in the initial stages of 
chemical synthesis in the primordial Universe chemistry 
(GP98, Lepp et al. 2002, Dalgarno 2005, Hirata & Pad- 
manhaban 2006). It represents one of the reagents for the 
process: 

HeH+-|-H^ H;^+He 

which, together with H+H"*" radiative association, is the 
main formation pathway for the molecular hydrogen cation. 
The equilibrium, permanent dipole moment of HeH^ is ap- 
proximately equal to 1.66 D (Pavanello et al. 2005). Engel 
et al. (2005) computed line lists for different HeH"*" isotopo- 
logues (^HeH+,''HeH+, ^HeD+,''HeD+), and data are avail- 
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Figure 2. HD"*" radiative cooling function: Circles: present cal- 
culation; red dashed curve: fit by Glover & Abel 2008; green solid 
curve: fit given in Table [Tl 
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Figure 3. HeH"*" radiative cooling function contributions. Circles: 
•^11011+; squares: ^HcD+; green solid curve: ''HeH+ fit (Table[TJ; 
blue dashed curve: 3HcD+ fit (Tablc[TJ. 

able electronicallj{3 both for ro-vibrational levels and Ein- 
stein coefficients. The cooling functions for all the isotopo- 
logues are shown in Figs. [3]and|4l 

2.4 LiH 

LiH represents the primordial molecule with the largest 
permanent equilibrium dipole moment (« 5.88 D, 
Hertzberg 1978). Several studies have modelled in detail 
lithium kinetics in the primordial Universe chemistry (e.g. 
Bougleux & Galli 1997, Lepp et al. 2002), also explicitly 
evaluating the cooling properties of lithium hydride due to 
coUisional processes. The cooling properties have been de- 
scribed by GP98, where a fit of the coUisional cooling func- 
tion in the low density limit was given. 

^ exomol - Molecular Line lists for Exoplanet Atmospheres, 
www.cxomol.com 
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Figure 4. HeH"*" radiative cooling function contributions. Circles: 
^HcH+; squares: ^HeD"*"; green solid curve: ^HeH"*" fit (Table[TJ; 
blue dashed curve: ''HcD+ fit (Tabic [TJ. 



''LiH has a singlet X ^E"*" ground state with an experi- 
mental equilibrium distance = 3.015 ao and a dissociation 
energy De = 20,287.7 cm"^ (Stwalley & Zemke 1993). 

In virtue of its small number of electrons LiH has 
been the object of many quantum mechanical studies, 
among which we may mention the pioneering one b y 



These studies, however, do not provide in a readily usable 
format the potential energy curve and dipole moment curve 
we need to compute the cooling function. We therefore in- 
dependently computed the required quantities. 



2.4-1 Computational detail 

The potential energy and dip ole mome nt curves have been 
calc ulated using the program ImolproI . supplemented with 
the ImrccI package. Descriptions of the quantum chemical 
methods used, along with referen ces to the original papers , 
can be found, e.g., in the review bv lLodi fc Tennvsonl(|2010l ). 

The basis sets used for hydrogen belon g to aug-cc- 
pVnZ, (n = T, Q, 5) family bv lDunninj ^98^)- For lithium 
we use d the recent aug - cc-pw CnZ(-DK) (n — T, Q, 5) basis 
sets bv lPrascher et ahl l|2010t ). Basis sets with the DK suf- 
fix were especially devised to be used with the relativistic 
DKH hamiltonian (see below). Computationally it would be 
entirely possible to use larger basis sets but, regrettably, no 
larger ones are available for lithium. 

The restricted Hartree-Fock (rhf) method applied to 
LiH wrongly dissociates to the ionic limit Li+(ls2)-FH~(ls2) 
instead of the correct neutral pathway Li(ls'^2s^)-|-H(ls^). 
This behaviour is a direct consequence of the rhf model, 
which upon dissociation restricts molecular fragments to be 
spin-singlets (see e.g. Lodi & Tennyson 2010). Because the 
dissociation limit is wrong, methods such as CCSd(t) which 
rely on the rhf solution being a good approximation to the 
exact wave function may be expected to experience prob- 
lems at long bond lengths. In fact, the CCSd(t) module in 
MOLPRO does not converge for r greater than about 6 ao 
(using the default settings). This may be slightly suprising 



because, as LiH has only 2 valence electron, method such 
as CCSd(t) are formally equivalent to full configuration in- 
teraction (fci) in the valence space and should therefore be 
quite accurate. It is the non-linearity of the equations to be 
solved, together with the poor starting guess, that result in 
these numerical problems. 

We calculated for six geometries all-electron full con- 
figuration interaction (fci) reference energies in the aug-cc- 
pwCQZ basis sets using the CCSDTQ code of mrcc. Com- 
putation of each FCI point involves about 36 x 10® configu- 
ration state functions and took between 30h and 50h on a 
4-core Intel Xeon 5160 workstation with 16 GiB of memory. 
We compared these FCI values with several single-reference 
and multi-reference methods; the results of the comparisons, 
not discussed here, indicate that a very good balance of 
speed and accuracy is provided by the internally-contracted, 
multi-reference averaged coupled pair functional (ic-ACPf) 
method based on a full-valence (2-electron, 5-orbital) com- 
plete active space (8 reference configurations). This method 
was therefore used for further calculations. With respect to 
the FCI energies the ic-ACPF method predicts a dissociation 
energy too low by about 6.0 cm~^. Computation of a IC- 
ACPF energy in the aug-cc-pCw5Z-DK basis set took about 
4 min on the same Xeon 5160 machine. 

To account for scalar-relativistic effects we considered 
the use of the traditional mass- velocity one-electron Darwin 
operator (mvdI) or, alternatively, of the fourth-order Dirac- 
KroU-Hess hamiltonian (dkh4). Relative energies computed 
by the two methods agree to better than ±0.05 cm~^. We 
chose the dkh4 approach. Relativistic corrections are small, 
affecting relative energies only by ±2 cm^^. The contribu- 
tion of relativistic corrections to the dissociation energy is 
—0.5 cm^^. Finally, we remark that in our tests differences in 
relative energies between DK-contracted and aug-cc-pwCnZ 
basis sets are negligible, less than ±0.05 cm~^ for energies 
up to dissociation. 

The energy and dipole curves were computed for 300 
uniformely-spaced grid points for r — 1.00 ao to r = 16.00 ao 
with spacing 0.05 ao. The final energy values were obtained 
for each grid point by basis-set extrapolation of the n =T, 
Q and 5 values; following common practice (see, e.g., Lodi 
& Tennyson 2010) the CAS-SCF energies were extrapolated 
using the functional form En — Eoo + Ae^"" while the dif- 
ferences between ic-ACPF and CAS-SCF energies were extrap- 
olated using E„ = _Eoo + A/n^. 

A minor complication arose because, for all basis sets 
employed, the CAS-SCF energies experience a discontinuity 
jump of about 6 cm^^ around r = 11 ao if RHF orbitals are 
used as starting guess; as a result the ic-ACPF energies also 
experience a discontinuity jump of about 1.2 cm~^. This 
problem is undoubtedly another consequence of the poor 
RHF starting guess; it was circumvented by starting calcu- 
lations at small bond lengths and then concatenating sub- 
sequent calculations for increasing r so that the CAS-SCF 
orbitals from the preceeding geometry were used as starting 
guess. 

The final potential energy curve has an equilibrium 
bond length Ve = 3.014 ao and a dissociation energy 
De — 20,294 cm^^; the dipole moment curve has an equi- 
librium dipole moment of 2.293 a.u., reaching a maximum 
of 3.003 a.u. for r — 5.180 ao. We may assign to De an 
uncertainty of the order of ±15 cm""^, mainly due to ba- 



Radiative cooling functions for primordial molecules 5 








, 1 1 1 , 1 , 1 , 1 , 1 1 [ 1 [ 1 [ 1 [ , 1 1 1 1 _ 






PEC 




-5000 










10000 










15000 




1 / \ DMC 

\ / V 






20000 







_ 1 , 1 , 1 


1 



Bond length in hnlirs 



modest dissociation energy De of 1130.5 cm"'^, an equilib- 
rium bond length = 4.130 ao and a permanent equilib- 
rium dipole (measured in the molecular centre-of-mass ref- 
erence frame) of 0.2773 a.u. (0.7048 D). 

Very specialised, high-accuracy methods a re applicable 
to_LiH+ due to its small size. Very recently iBubin et al.l 
presented non-adiabatic, relativistic-corrected values 
for the first five rotationless ( J = 0) vibrational states with a 
stated accuracy of better than 0. 1 cm~ ^ . Among other recent 
studies tre ating LiH ^ we m ay mention iGadea fc Leiningerl 
(|2006h and iMagnied { 2004h . However, again as in the case 
of LiH, because of the different focus these recent studies 
do not provide in a usable-format both the potential energy 
curve and the dipole moment curve we need to compute the 
cooling function. We therefore independently computed the 
necessary quantities. 



Figure 5. LiH potential energy curve (PEC) and dipole moment 
curve (dmc). The two curves use different scales for the y axis; 
the left-hand side legend refers to the PEC, the right-hand side 
one to the DMC. 
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Figure 6. LiH radiative cooling function. Circles; present calcu- 
lation; green solid curve: fit given in Table [T] 

sis set incompleteness. For c omparison, the non-r e lativis tic 
'exact BO' value quoted by ICooper fc Dickinson! (|2009t ) is 
De = 20, 298.8 cm~^. The PEC and the DMC are plotted in 
Fig.H 

2.5 LiH+ 

LiH+ represents the most abundant molecular species con- 
taining lithium at low redshift. Indeed, while in the very be- 
ginning phases LiH is much more abundant than its cation, 
this behaviour reverses at « ~ 20, as predicted by Dalgarno 
& Lepp (1987) and subsequently found by kinetics models 
(e.g. GP98, Bougleaux & Galh 1997). The reason of this 
behaviour is the presence of a residual ionization fraction 
which enhances one of the main formation channel, namely 
ion-atom radiative association. 

So far LiH"*" has not been characterised experimentally. 
The present calculations, discussed below, predict a very 



2.5.1 Computational details 

Like the LiH case discussed in Section r2.4.1l we used the aug- 
cc-pwCnZ-DK {n — T, Q, 5) basis sets and we accounted 
for scalar relativistic effects by speci fying the dkh4 hamil- 
tonian. Calculations were done using ImolproI . 

The restricted open-shell Hartree-Fock method (rohf) 
applied to LiH""" correctly dissociates to Li^ -I- H. Single- 
reference methods are therefore expected to fare very well, 
in contrast to the LiH case discussed in Section [2.4. II 

As a preliminary test we compared using the aug- 
cc-pwCTZ-DK basis set 38 full configuration interac- 
tion (fci) relative energies in the range 3.00-11.00 ao 
with energies computed by various other methods. 
Here is a brief summary of the comparison, in the 
form METHOD/error; rohf/98.9 cm~^; CISd/1.8 cm"^; 
UCCSd/1.7 cm"^; CISD-|-p/1.4 cm~^; CISD+q/0.6 cm"^; 
UCCSd(t)/0.3 cm^^. The reported errors are one half of the 
non-parallelity error, which is defined as the difference be- 
tween the maximum and the minimum of the modulus of the 
energy differences between two methods (Li & Paldus 1995). 
As expected, uccsd(t) produces extremely accurate results; 
also note that this method is asymptotic to FCI upon dis- 
sociation as it is size-extensive and exact for two-electron 
sy stems. W e therefore used rhf-uccsd(t) as implemented 
in ImolproI for further calculations. 

Because LiH""" is a charged species the value of the 
dipole moment depen ds on the choice of the origin of the 
axes. As discussed bv lBunker fc JensenI ()l998h . the correct 
choice of the origin is the molecular centre-of-mass. mol- 
PRO by default does set the origin at the centre-of-mass but 
uses isotopically-averaged nuclear masses, which is not ap- 
propriate in our case. We therefore specified the following 
masses: m(Li+) = 7.015455 u and m(H) = 1.0078250321 u. 
These are atomic masses, to partially account for the elec- 
tron contribution to the centre-of-mass coordinates; m(Li"^) 
was obtained by subtracting one electron mass to the atomic 
mass of Li. 

With this choice of masses and taking the molecule 
to lie along the x axis the coordinates of the nuclei for 
an inter-nuclear distance r are: 2;(Li^) = 0.125613r and 
2;(H) = -0.874387r. The dipole moment M{r) is therefore 
asymtotic to M{r) — > 0.125613r in this reference system. 

We present as figure [7] a plot of the PEC and of the 
quantity 0.125613r — M(r), namely the difference between 
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Figure 7. LiH"*" potential energy curve (PEC) and the quantity 
0.125613r — M{r) (see text for details). The two curves use dif- 
ferent scales for the y axis; the left-hand side legend refers to the 
PEC, the right-hand side one to the DMC. 
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Figure 8. LiH"*" radiative cooling function: circles: present cal- 
culation; green solid line: fit given in Table [T] 

the asymptotic dipole and the actual one. We chose not to 
plot directly M{r) as, because of the asymptotic form, it 
would very nearly appear as a straight line. 

We report here below the rotationless (J = 0) vi- 
brational terms Eiy -I- 1) — E{v) toget her with the differ- 
ence with the very accurate values by iBubin et al.] (|201ll ) 
for V = 0,1,2,3,4: 354.59(0.44); 261.79(0.78); 170.00(0.59); 
89.82(-0.07); 35.26(0.18). Our computed spectroscopic val- 
ues are therefore very accurate, with errors of less than 
1 cm-^ 
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Figure 9. H2D"'" radiative cooling function: circles: calculation 
by Sochi &: Tennyson 2010; green solid line: fit given in Table [T] 

I Glover fc Savi 3 (|2009l ). The former work gives a fit for the 
radiative cooling function as a function of temperature based 
on the ab initio line list of Neale et al. (1996). 

2.7 H2D+ 

In the deuterium chemistry of primordial Universe, H2D"'" 
represents the most complex t riatomic molecule usua lly in- 
troduced (jStancil et al.lll998l : IVonlanthen et al.l 120091 ') . Re- 
cently, Sochi & Tennyson (2010) have calculated a compre- 
hensive line list of H2D"^ frequencies and transition proba- 
bilities. Table [T] provides an improved fit for the radiative 
cooling function at low temperatures. 



3 CONCLUSIONS 

In the present work, the cooling functions of the most abun- 
dant molecular species in the primordial Universe have been 
studied; calculations of the radiative contributions have been 
taken into account, under the hypothesis of local thermal 
equilibrium distributions of ro-vibrational levels. Radiative 
cooling functions are shown for each molecule. Analytic fits 
to each of these functions were obtained in the form: 

JV 

log.o W = ^a„(logior)". (7) 

n = 

Coefficients for these fits, a-n, are provided in Table [1] The 
fits are valid only in the temperature range specified in the 
figures. Potential energy and dipole moment curves calcu- 
lated in the present work, together with the computed tran- 
sition wavenumbers and Einstein coefficients, can be down- 
loaded from www.exomol.com. 



2.6 H+ 

is one of the most interesting molecular ion in 
astrophysics, especially for its role in the interstellar 
medium chemistry. Its cooling properties have been re- 
ce ntly considered, bot h in planetary atmosphere conditions 
bv lMiUer etahl (|2010l ) and in the primordial Universe by 
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Table 1. Radiative cooling function fits. 



Molecule Coefficients Molecule Coefficients Molecule Coefficients 



HD 


N= 


=5 


ao 


= 


-55.5725 


HD+ 


N= 


7 


ao 


= 


-6.04917 


3HeH+ 


N 


=7 


OQ 




-10.2807 








0.1 


= 


56.649 








ai 


= 


-60.0312 








Ol 


= 


-62.9415 








12 


= 


-37.9102 








12 


= 


98.8361 








02 


= 


118.348 








a-i 




12.698 








0-3 




—77.5575 








03 




— 97.3721 








a4 




-2.02424 








04 


— 


33.4951 








a4 


— 


42.5517 








0.5 


= 


0.122393 








15 


= 


-8.07092 








05 


= 


-10.1946 




















ae 


= 


1.01514 








06 


= 


1.26335 




















07 


= 


-0.0519287 








07 


= 


-0.0633433 


*HeH+ 


N= 


=7 


oo 


= 


-7.58736 


3HeD+ 


N= 


7 


ao 


= 


26.2717 


''HeD+ 


N 


=7 


00 


= 


28.5384 








Ol 


= 


-68.2966 








0.1 


= 


-168.493 








Ol 


= 


-172.458 








02 


= 


122.847 








a2 


= 


247.288 








02 


= 


249.811 








03 




on AAA 

— yy .444 








13 




— 164. zyy 








03 




— io4. / 00 








04 




43.1409 








a4 


= 


77.0755 








04 


= 


76.9479 








05 




-10.3034 








05 




— 18.2012 








05 




— 18.1306 








06 




1.27565 








ae 




2.26293 








06 




2.25209 








07 




-0.0639846 








17 




-0.11512 








07 




-0.114559 


LiH 


N= 


=6 


00 




-31.894 


LiH+ 


N= 


7 


ao 




-23.5448 


H2D+ 


N- 


=6 


OQ 




33.8462 








01 




34.3512 








ai 




9.93105 








01 




-188.249 








02 




-31.0805 








0.2 




-8.6467 








02 




253.037 








03 




14.9459 








03 




2.13166 








03 




-168.02 








04 




-3.72318 








04 




2.43072 








04 




59.3597 








05 




0.455555 








05 




-1.69457 








05 




-10.6334 








Og 




-0.0216129 








ae 




0.384871 








Og 




0.759029 




















07 




-0.0300114 




















fit bv Miller et al. C2010) 
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